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ABSTRACT 

Clumping in hot star winds can significantly affect estimates of mass-loss rates, the inferred evolution 
of the star and the environmental impact of the wind. A hydrodynamical simulation of a colliding 
winds binary (CWB) with clumpy winds reveals that the clumps are rapidly destroyed after passing 
through the confining shocks of the wind-wind collision region (WCR) for reasonable parameters of 
the clumps if the flow in the WCR is adiabatic. Despite large density and temperature fluctuations 
in the post-shock gas, the overall effect of the interaction is to smooth the existing structure in the 
winds. Averaged over the entire interaction region, the resulting X-ray emission is very similar to 
that from the collision of smooth winds. The insensitivity of the X-ray emission to clumping suggests 
it is an excellent diagnostic of the stellar mass-loss rates (M) in wide CWBs, and may prove to be 
a useful addition to existing techniques for deriving M, many of which are extremely sensitive to 
clumping. Clumpy winds also have implications for a variety of phenomena at the WCR: particle 
acceleration may occur throughout the WCR due to supersonic MHD turbulence, re-acceleration at 
multiple shocks, and re-connection; a statistical description of the properties of the WCR may be 
required for studies of non-equilibrium ionization and the rate of electron heating; and the physical 
mixing of the two winds will be enhanced, as seems necessary to trigger dust formation. 

Subject headings: hydrodynamics — stars: individual: WR 140 — stars: mass-loss — stars: winds, 
outflows — stars: Wolf-Rayet — Xrays: stars 



1. INTRODUCTION 

There is now considerable observational evidence for 
a high degree o f structure, or "clumping", in hot ste l- 
lar winds (e.g. iMoffat et all 119881 : iLepine et all 12000( 1 . 
Clumping affects diagnostics of mass-loss rates which are 
sensitive to the square of the density (e.g., free- free radio, 
infra-red continuum, Ha). If not taken into account, the 
inferred mass-loss rates may be substantially above their 
actual values. For Wolf-Rayet (WR) stars, the assump- 
tion of a smooth wind typically leads to a factor of 3 
overestimate. Recent studies are now indicating that the 
winds of main sequence stars may be even more struc- 
tured, and downward revisions in M by factors of 3 — 10 
or more have been suggested (e.g.. IBouret et all 120051 : 
IFullerton et~aTI 120061 : IPuls et all |2006j). Such substan- 
tial reductions have a dramatic effect on the evolution 
and environmental impact of massive stars, and may also 
be needed to explain th e near-symmetry of X-ray lines 
(|Owocki fc Cohenll2006Tl . 

Several methods for determining mass-loss rates which 
are not sensitive to p 2 are unique to binary systems. 
First, if the period of the binary is short enough, and 
the mass-loss rates high enough, M can be determined 
from the observed change in the orbital period. Though 
this has been applied to the WR binary V444 Cyg (e.g. 
lAntokhin et al.lll995f ). it has limited applicability. Ob- 
served changes in polarization and atmospheric contin- 
uum eclipses can also be used to determine mass -loss 
rates (|St. -Louis et al.l 1988: La montagne et al.lll996( l. 

Measurements of M have also been made by com- 
paring the observed X-ray flux arising from the WCR 
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to predict ions from hydrodynamical models of this in- 
teraction ([Stevens et ai1ll996l : lPittard fc Corcoran1l2002l : 
Pittard & Dougherty 2006). In all of these models homo- 
geneous winds were assumed, but this work shows that 
structured winds are rapidly smoothed out in the WCR 
of wide binaries: hence, the X-ray emission, despite be- 
ing sensitive to p 2 , is an excellent diagnostic of the actual 
mass-loss rates. Further implications of wind structure 
on the interaction region are discussed in § [4] 

2. CLUMP- WCR INTERACTION 

The survival time of a clump within the WCR depends 
predominantly on its size and density contrast to the 
mean flow. The low amplitude of detected variability ar- 
gues for a large number of clumps, each of which is likely 
to have a quite small spatial scale (e.g. lEversberg et all 
1998; Mar chenko et alll2006l and references therein). If 
the interclump medium is devoid of material, the den- 
sity contrast of the clumps relative to the corresponding 
smooth flow is inversely related to the volume filling fac- 
tor, f v , which the clumps occupy. f v initially decreases 
with radius (i.e. the wind becomes increasingly clumpy) , 
reaches a minimum at around 10 — 15 R*, and then de - 
clines as the wind slowl y smooths out (IPuls et al.|[ 2006l 
At very large distances, iRunacres fc Owockil (|2005h find 
that / v ~ 0.25. 

Clumps which pass into the WCR lose mass primarily 
through dynamical instabilities (e.g. Kelvin-Helmholtz). 
The destruction timescale of no n-radiative clump s can be 
parameterized as t^ — er c /v s ([Klein et al.lll994h . where 
r c is the clump radius and v s is the shock velocity (for 
the stationary shocks in CWBs, v s is equal to the normal 
component of the pre-shock wind speed), e « 3.5 for a 
density ratio between the clump and inter-clump wind of 
order 10 — 100. A detailed review of clump destruction 
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processes, including the effects of radiative co oling and 
magnetic fields, can be found in[Pittard| (j2006|). 

Consider a clump moving along the line of centers be- 
tween the stars. As the clump passes through one of the 
shocks confining the WCR it is decelerated less than the 
inter-clump material, and pushes the confining shock into 
the WCR. A lower limit on the timescale for the clump 
to half cross the WCR is t cd = Ar/v s , where Ar is half 
the distance between the confining shocks. Clumps will 
be destroyed before they reach this half-way point when 
td/tcd = er c / Ar < 1. Since it is expected that r c < Ar, 
clumps should be rapidly destroyed within the WCR. 

3. HYDRO DYNAMICAL SIMULATIONS 

The long-period CWB WR 140 forms the basis of 
this investigation into the effects of clumpy winds 
on the interaction region. WR 140 is the achetype 
of long-period CWBs, exh i biting dramatic variations 
in its X-ray (Poll ock et all 120051 ) and radio emission 
(|Doughertv et al.ll2005t ) modulated by its highly eccen- 
tric orbit. The emission fr om radio to TeV energies ha s 
recently been modelled bv lPittard fe Dougherty! (|2006). 

There have been two previo us studies of the c lumpy 
colliding winds in WR140. IWalder fc Folinil (|2002f ) 
examined a simulation with the stars near perias- 
tron. WR material both in and between the clumps 
rapidly cools, and is c o mpres sed to high densities. 
lAleksandrova fc Bvchkovi ((2000:) considered the transi- 
tion between periastron and apastron in an attempt to 
explain the variation of the X-ray flux with the stellar 
separation. 

In this work the effect of clumpy winds on the WCR 
is examined when the stars are at apastron (i.e. the 
stellar separation D — 4.7 x 10 14 cm). The mass-loss 
rates and terminal wind speeds adopted are: AfwR = 
4.3 x 10~ 5 Moyr" 1 , M = 8.0 x 10" 7 Mgyr" 1 , 
Voo wr = 2860 kms™ 1 , v^o — 3100 kms -1 . In order to 
focus attention on the apex of the WCR, which is bent 
sharply around the O star, the WR star is positioned off 
the 2D hydrodynamic grid, which is axisymmetric and 
contains 1650 x 660 cells. The clumps are added to the 
flow in annuli around each star at specific time intervals. 
For simplicity the clumps are assigned the wind terminal 
speed, and have radii proportional to their distance from 
their star. The clump radius, r c , is 0.005 r in the WR 
wind, and 0.02 r in the O wind, so that the clumps have 
> 10 cells across their radius when interacting with the 
WCR. The clumps are given a density contrast of 10 with 
respect to the interclump medium, and the clump and 
interclump medium are assumed to contain equal mass 
(hence f v = 1/11). The inter-clump medium is perfectly 
smooth. A simulation with smooth winds reveals that 
the width of the WCR on the line between the centers 
of the stars is 0.0695 D. Hence t d /t cd w 0.5 (0.25) for 
the clumps in the WR (O) wind, and the WCR should 
quickly smooth out these inhomogeneities. 

The hydrodyn amical code is second-order accurate in 
space and time (|Falle fc Komissarovl 119960 . Optically 
thin radiative cooling is included, though radiation losses 
are negligible. Heat conduction is not explicitly included, 
and is likely to be strongly inhibited by the magnetic 
field revealed by the synchrotron emission. The clumps 
are initially spherical, but the velocity dispersion in the 
radial and transverse directions may differ. vg/v r w 4 for 




3x10' 4 



3.5x10 14 



4x10 14 4.5x10 14 
Distance (cm) 



Fig. 1. — a) Logarithmic density plot (logio gem -3 ) from an 
axisymmetric simulation of colliding clumped winds in WR 140. 
b) As a) but showing the Mach number of the gas. Supersonic gas 
is white. Due to the pressure gradients within the WCR, the gas, 
on average, becomes supersonic as it leaves the system. 



the clumps in the O wind, so that they naturally "pan- 
cake" as they are advected by the wind. Recent models 
of X-ray wind absorption s upport radially comp ressed 
rather than spherical clumps (jOskinova et alJ feOOG) . The 
density contrast of the clumps is also reduced by numer- 
ical diffusion. These two effects reduce the survival time 
of the clumps, but are not thought to be significant. 

A density plot of the WCR is shown in Fig. QJi). It is 
immediately evident that the process of clump destruc- 
tion induces a multitude of large and small scale motions 
within the WCR. This "turbulence", some of which is 
supersonic (Figs. [Ib and[2H), puffs-up the WCR, so that 
its volume is 25% larger than the homogeneous case. 
Temperatures higher than in the smooth-wind case occur 
when the WCR expands upstream into the relatively low 
ram-pressure of the interclump medium, and also behind 
the bow-shocks driven ahead of the clumps. In contrast, 
the interiors of clumps are initially heated to significantly 
lower temperatures as shocks driven into them are much 
slower (shocks driven into clumps in the WR wind near 
the line of centers have speeds ~ 1000 kms , and pro- 
duce temperatures of ~ 4x 10 7 K). Thus, there is a wider 
distribution of temperatures within the WCR than in the 
smooth winds case (Fig. ^>). 

As cool material within recently-shocked clumps moves 
deeper into the WCR it is heated by secondary shocks 
and through mixing with hotter plasma. At any partic- 
ular instant, the majority of the mass within the WCR 
is heated to temperatures similar to those that exist in 
the smooth winds simulation: the mass-weighted mean 
temperature of gas within the WCR (and on the grid) is 
2.4 x 10 8 K and 2.0 x 10 s K for the smooth and struc- 
tured winds cases respectively. The average density in 
the WCR is similar to the smooth winds case, and is not 
markedly different with 90% of the mass in clumps. 

Although the WCR is clearly not smooth, the wind- 
wind collision decreases the small-scale structure in the 
winds. Clump disruption is mostly through the devel- 
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Fig. 2. — a) Comparison of mass as a function of Mach num- 
ber from simulations with homogeneous (solid line) and clumpy 
(dashed line) winds, b) As a) but as a function of temperature. 

opment of large scale perturbations, while the mixing of 
the clump and interclump material is affected by smaller 
scale motions. The global nature of the simulation in- 
evitably leads to relatively poor spatial resolution at the 
scale of individual clumps, and is likely to enhance the 
mixing due to non-negligible numerical diffusion. In con- 
trast, the lack of a third dimension in these simulations 
(the imposition of axisymmetry means that the clumps 
are donut shaped) should slow the rate of clump destruc- 
tion as there is one less degree of freedom for dynamical 
instabilities. Resolution tests and comparison with previ- 
ous work on shock-cloud interactions indicate that these 
effects do not have a serious impact on the results. 

In reality, clumps will possess a variety of density con- 
trasts and sizes - large, dense clumps will survive longer 
as distinct entities within the WCR than smaller, less 
dense clumps, though d enser clumps may on average be 
smaller and vice- versa (Moffat 1994). If id/ied >> 1, 
clumps could, in theory, pass completely through the 
WCR and into the pre-shock wind of the companion star, 
but this is unlikely to occur in the wide, adiabatic sys- 
tems considered here. 

3.1. Determining Mass-loss Rates 

Since the wind structure is rapidly smoothed, it is not 
surprising that the X-ray emission is similar to the ho- 
mogeneous case (Fig. [3]). The mean continuum flux at 
1 keV is only 17% higher in the clumpy winds simulation, 
which translates into an 8% overestimate of the mass-loss 

. 2 

rates (/ x oc M ). Clearly, there is potential to use the X- 
ray emission from the WCR as a diagnostic of the stellar 
mass-loss rates. 

However, the structure within the WCR and the re- 
sulting X-ray emission depends on the properties of the 
clumps. Further simulations (Pittard, in preparation) 
reveal that the emission measure (J n 2 alV) appears to 
increase with the clump density contrast: an identical 
simulation with a density contrast of 100 yields a 56% 
overestimate of M for each star, but the density contrast 
is unlikely to be this high in wide CWB systems. In- 
stead, the major sources of uncertainty are likely to be 
the abundances and the wind momentum ratio, 77 (as the 
emission is nearly degenera te between 77 and the M's - see 
iPittard fc Dou gherty 2006). Nevertheless, X-ray derived 
mass-loss rates should be accurate to within a factor of 
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Fig. 3. — Comparison of the X-ray emission from simulations 
with homogeneous (thick line) and clumpy (thin lines) winds. The 
latter is slightly softer, and shows only weak variability (< 5%), 
despite the fact that the hydrodynamic grid is not large enough 
to capture all of the emission (note the lack of lines b elow 1 keV). 
The emission was calculated using the MEKAL code l|Mewe et al.l 
[1993, and references therein), for optically thin thermal plasma, 
ionization equilibrium, and identical ion and electron temperatues. 

2. Substantially higher precision may be obtained if 77 
and the wind abundances are strongly constrained. 

The level of X-ray absorption by the intervening 
wind(s) can also, in principle, be used to determine mass- 
loss rates. For the wide binaries considered here the 
clumps are likely to be optically thin, so that the op- 
tical depth will be identical to the smooth winds case. 
However, again there is a degeneracy with 77 if the 
line-of-sight into the system is not well constrained, in 
which case it may be best to simply match the X-ray 
flux at energies above thos e susceptible to absorption 
(jPittard fe Doughertvll2006T) . 

4. FURTHER IMPLICATIONS 

The interaction of clumpy, as opposed to homogeneous, 
winds has implications for a variety of phenomena which 
occur at or within the WCR. 

4.1. Particle Acceleration and Synchrotron Emission 

The radio synchrotron emission from wide CWBs has 
recently been modelled assuming that relativistic elec- 
trons are created by diffusive sh ock acceleration (PSA) 
at the shocks confining the WCR (iDoughertv et al1l2003l : 
IPittard et al.ll20q6t IPittard fc Dougherty! 120061) . In the 
case of WR 140. IPittard fc Dougherty! (|2006f) find that 
the non-thermal (NT) electron energy distribution is 
harder than the canonical DSA value, i.e. p < 2. While 
there are many possible explanations, the results pre- 
sented above hint at several mechanisms. First, NT 
particles may be accelerated via the second-order Fermi 
process resulting fro m clump-induced MHD t urbulence 
within the WCR (e.g. lScott fc Chevalieilfl975T) . Second, 
particles accelerated at the confining shocks may be re- 
accelerated at m ultiple (weak) shocks within the WCR 
(|Schneiderill993h . A third possibility is magnetic recon- 
nection, which probably occurs throughout the volume of 
the turbulent WCR, and not just at a hypothetical con- 
tact discontinuity. Reconnection may also provide addi- 
tional energy for generating and maintaining the mag- 
netic fluctuations which drive stochastic acceleration. If 
any of these mechanisms is dominant, models which im- 
pose particle acceleration only at the confining shocks 
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will need to be revised accordingly. In addition, the 
magnetic field within the WCR is likely to be highly 
"tangled" due to the turbulent motions of the gas in- 
side. This was a central assumption in the recent models 
mentioned above, and means that the asymmetric radio 
lightcurve of WR 140 cannot be explained by the angular 
dependence of the synchrotron emission process. 

4.2. Plasma Timescales 

There are a variety of timescales within the WCR 
whic h may be modified if t he stellar winds are clumpy 
(e.g. JWalder fc Folinill2003) . In wide CWB systems, the 
shocks are collisionless, and the electron and ion temper- 
atures may differ (an effect which appears to be sensi- 
tive to the shock speed, e.g. iHwang et al . 2002). Equili- 
bration subsequently occurs through Coulomb collisions, 
proceeding faster where the density is higher, such as in 
material within the shocked clumps. The spectral hard- 
ness is sensitive to this process, but the overall flux is 
not dramatically changed, and it is unlikely to have a 
significant impact on mass-loss rate determinations. 

Ionization equilibrium may also take a significant time 
to occur. Direct evid ence for non-equilib rium ioniza- 
tion is seen in WR140 (Pol lock et al.1120051 ). and models 
which assume ionization e quilibrium fail to rep roduce ob- 
served X-ray line profiles (|Henlev et al.l 12005). Material 
originally within clumps will ionize quicker than inter- 
clump material, but to lower stages due to the reduced 
post-shock temperatures. The relevant timescale to ob- 
tain the highest ionization stages then becomes the speed 
at which further heating occurs in the downstream flow, 
through compressions, shocks, and mixing. The details 
are again sensitive to the clump properties, but the con- 
tinuum emission, and thus estimates of M, should not 
be strong affected. 

4.3. Dust Formation 

The highly turbulent interior of the WCR shown in 
Fig. [1] enhances the mixing between the two winds. Such 
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mixing may be necessary in order that carbon-rich WR 
material and hydrogen-ric h O-star material can form 
dust within the WCR (e.g. lWalder k Folinil[200l . 

5. SUMMARY 

The interaction of clumpy stellar winds in massive bi- 
nary systems creates a highly turbulent wind- wind colli- 
sion region in the adiabatic limit. The lifetime of clumps 
within the WCR depends on their density contrast and 
size. Clumps with a density contrast of 10 and radii a 
few times smaller than the half- width of the WCR on the 
line of centers between the stars are rapidly destroyed. 
Material originally within the clumps is then vigourously 
mixed into the surrounding flow. 

The stochastic impact of clumps on the WCR and their 
subsequent destruction creates significant density and 
temperature fluctuations within the WCR, but the global 
X-ray emission can be remarkably similar to the smooth 
winds case. The X-ray emission is then an effective, 
clumping- independent, measure of the stellar mass-loss 
rates, especially if the wind momentum ratio is known. 
The small number of CWB systems suitable for such an 
analysis is countered by the potential accuracy that can 
be obtained. Each time this method has been used to 
determine mass-loss rates, values lower than those in the 
literature were inferred. Although only a small part of 
parameter space is explored here, any reasonable param- 
eters for the clumps in the wide, adiabatic, systems con- 
sidered here should lead to their rapid dissolution. 

Turbulence and weak shocks within the WCR provide 
mechanisms for obtaining hard NT particle spectra. The 
timescale to obtain high ionization stages may be con- 
trolled by the rapidity of clump destruction, and en- 
hanced mixing of the winds will aid dust formation. 
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